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^ ■ Abstract 

3^ I We complete the construction of raising and lowering operators, given in a 

previous work, for the orthogonal polynomials of hypergeometric type on non- 
homogeneous lattice, and extend these operators to the generalized orthogonal poly- 
nomials, namely, those difference of orthogonal polynomials that satisfy a similar 
^ ' difference equation of hypergeometric type. 
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Recently we have presented a paper on the raising and lowering operators for the orthogo- 
Ph' nal polynomials (OP) of hypergeometric type [5] (in connection with the factorization 
r-| ■ method defined by Hull and Infeld). That paper covered only OP of continuous and 
^ ! discrete variable on an uniform lattice, as well as orthonormal functions of continuous 
^ I and discrete variable. 

■ In this work we continue the construction of raising and lowering operators for OP 

k>( ' non-homogeneous lattice. The starting point is also the Rodriguez formula and the 

. fundamental properties of OP given by Nikiforov and collaborators [1], [18] that include 
5^ ! the q-analog of classical OP of discrete variable. 

Atakishiyev and collaborators [13] [14] extended the classification of Nikiforov to OP of 
discrete variable defined by Andrews and Askey and proved that they satisfy a difference 
equation only in the case of x{s) linear, quadratic q-linear and q-quadratic. 

The construction of raising and lowering operators on non-uniform lattice was also 
worked out by Alvarez- No darse and Costas-Santos [15] [16] for the lattice x(s) = Cig** + 
C2q~'' + C3. 

Similar work was carried out by Smirnov for OP of hypergeometric type on homo- 
geneous lattice [20] and non- homogeneous lattice [21] although the raising and lowering 
operators are defined with respect to two indices: n, the order of polinomials and m, 
the order of difference derivative of polynomials. In our work we define the raising and 
lowering operator with respect to one index only, n or m. 
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In order to complete the classification of the OP of hypergeometric type we include the 
generalized classical OP that satisfy also a difference/differential equation of hypergeo- 
metric type. 

Since all classical OP of discrete variable lead in the limit to the corresponding OP of 
continuous variable, we start in section 2 with the raising and lowering operators for gener- 
alized OP of continuous variable with respect to the index n, using Rodrigucs formula. In 
section 3 we repeat the same construction for generalized classical OP of discrete variable 
on hmogeneous lattice. In section 4 we extend the construction to classical OP of discrete 
variable on non-homogeneous lattice in the general case, when x{s) = cis^ -|- C2S + C3 or 
x{s) = ciq^ + C2q~^ + C3. 

In section 5 we complete the picture with the construction of raising and lowering ope- 
rators for generalized classical OP on non-homogeneous lattice, that include the q-analog 
of classical OP of discrete variable. In all these cases, the raising and lowering operators 
are given with respect to one index, say, n, but the same operator can be considered, 
written in appropiate form, the raising and lowering operator with respect to index m. 

2 Raising and lowering operators for generalized 
classical orthogonal polynomials of continuous 
variable 

Let yn{x) be an orthogonal polynomials of continuous variable satisfying the differential 
equation [1] 



where (t{x) and t{x) are polynomials of at most second and first degree, respectively, and 



It can be proved that the derivatives of yn{x), namely, yn \x) = Vmn{x),m = 0, 1, . . . 
n — 1, satisfy a similar equation: 



with Tm — t{x) + ma'{x) and jimn — —{n — m) (t' + "''^^'^ a") , m = 0, 1, . . . , n — 1 

We call these polynomials generalized orthogonal polynomials of hypergeometric type, 
some particular examples of which are the Legendre and Laguerre, Hermite, Jacobi gen- 
eralized orthogonal polynomials. [2] 

The polynomials of hypergeometric type satisfy an orthognality property with respect 
to the weight function p{x) 




(1) 




(2) 




(3) 



b 




(4) 



a 
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Similarly the generalized orthogonal polynomials satisfy 



/ 



where and are normalization constants. 
It can be proved [3] 

(n-l \ ~^ /n-1 

n ) ' = n ^fc^^ 

k=m J \k=0 

from which can be eliminated, therefore 

m—l 



A;=0 

where dgn ^^<i given in the tables [6] . 

The generalized orthogonal polynomials of hypergeometric type can be calculated from 
the weight function pm{x) = a{x)"^p{x), with the help of the Rodrigues formula: 

y (x) = "'^ " - in (x)\ (7) 

where 

m— 1 I m—l / \ \ 

^n.n=(-irn/^^"=r^n -^)- (8) 

11 (^_^)!11\^ ^ + 

The leading coefficients of the orthogonal polynomial yn{x) — anX^ + hnX^~^ + . . . can 
be calculated [4] 



n—l / , 



fc=0 ^ ^ 

hence, it follows, AnnBn = n!a„ 

We address ourselves to the construction of the raising and lowering operators for the 
generalized orthogonal polynomials using the Rodrigues formula as we did in a recent 
work [5] 

We have from (7) 



4 R ^n+1— m A ID jn—m 



^m,n+l -Bn+1 J _ ^^n^ C?" ™ ^^^^^ ^ _^_/ ^" 

cr™(a;)p(a;) 

^— <^ — Tn[x)vmn{x) + (n - m)- T^a[x)v^^[x, 

Bn+i ( n+1 n Xn+m .X .X + 1 n 



(10) 
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The right hand can be considered the raising operator that, when apphed to Vmn{x) 
gives a new polynomial of higher order 

In order to evaluate the lowering operator we need a recurrence relation for the gen- 
eralized polynomials. We write 



n+l 



fe=0 

b 

Ckn = / 'Vmk{x) X VmnPm{x)dx. (11) 



Prom the orthogonality condition (5) we deduce 



b 



Vmn{x) x^' Pm {x)dx = for T < u — m. 



a 



Since xp)^^ {x) is a polynomial of order k + 1 — m it follows that 
Ckn = if fc+l — m<n — m, or /c + l<n. Hence 

{x) + PnVmnix) + (x), (12) 

where O;^ C^+l^n) Pn ~ ^nm,) ^n—l,n 

The coefficients q;„, and 7„ can be expressed in terms of the squared norm and 
the leading coefficients and hn in yn{x). 

From (11) it can be proved that d^/^Ckn = d'^n^nk- 

Since cin-i = c„,„_i, 7n = c„_i_„, if we put /c = n — 1 we obtain 

72 ^2 

Cn-l,n"rn,n-l ~ '^n,n-l"rnn' 

hence 

~ m,n 

In Oin—1 



Introducing the expansion yn{x) — anX^ + ^ + . . . in (12) and comparing the 
coefficients of the highest terms, we have 

a„(n - m + 1) = a„a„+i(n + 1), 

hn{n - m) = CKn^n+in + /9„a„n. 



Hence 



an n-m + 1 

«n = — ; — , 13) 

Un+l n+l 

~ _ hn {n - m) hn+i n+l - m 

a„ n a„+i n + l 

a„_i n-m dl^ 

= ^ ■ ^^^^ 
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Substituting (9) in q:„ we obtain 



—Bn n — m + 1 Xn 2n 2n + 1 

«n = -B — j T 7 ■ (16) 



Hence (10) can be written 

an ^ Vm,n-rx{x) = I h±IL ZZiM Vmn{x) - Cr{x) v'^^{x)\ . (17) 

Inserting (12) in (17) we obtain 

In ^ Vrn,„-l{x) = |-^^^ + ^ ("^ ~ } ^"^"^"^^ + "^^"^^ ^"^"^"^^ ^^^^ 

The right hand side of (17) and (18) can be considered the raising and lowering oper- 
ators for the generahzed classical orthogonal polynomials with respect to the index n. 
All the constants $n,ln, ^n, ^'n can be calculated from the tables [6]. 
Now we define the orthonormalized function 

i^mnix) = d^:;^^^J Pm{x)Vmn{x), (19) 

hence 

(20) 

Multiplying (17) by d^^sj pm{x) and substituting (20) in (17) we get 

ttn^%^V'm,n+l(a;) = ^^y^^^^^V'mn(a;) + \rm-l{x)^mn{x) - (J {x)%l)'^^{x) = 



Similarly 

~ -^271 dm,n-l 



In 



In djjin 



^ L'^{x,n)ilJm,n{^) (21) 

'4^m,n-l{x) = |~~qr^"~r^ ~ '^") ~ ^'^rn-lix)^ 1pmn{x) + 

+ (^i^Wmni^) = L-(x, n)V'm,n(a:;) (22) 



that can be considered the raising and lowering operators for the generalized orthonormal 
functions ip^nix)- These operators are mutually adjoint with respect to the scalar product 
of unit weight. 

Following the same procedure as in [5] we can factorize the raising and lowering oper- 
ators as follows: 

L~{x, n + l)L~^{x, n) = ji{n) — (t{x)H{x, n) 
L'^{x, n)L~{x, n + 1) = p{n) — a{x)H{x, n + 1) 

where 

/ X X2n A2n+2 ~ ~ 

'^^ ' 2n 2n + 2 

and H{x, n) is the differential operator derived from the left hand side of (3) after substi- 
tuting 4>mn{x) instead of Vmn{x). 

Notice that the factorization of the raising and lowering operators is defined in a basis 
independent manner, which is equivalente to the Infeld-HuU method. 
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3 Raising and lowering operators for generalized 
classical orthogonal polynomials of discrete 
variable on uniform lattice 

Let Unix) be an orthogonal polynomial of discrete variable satisfying the difference equa- 
tion [7] 

a{x) AV Vnix) + t{x) Aynix) + A„y„(a;) = 0, (23) 
where a{x) and t{x) are polynomials of at most of second and first degree, respectively, 

Xn = -n(T' + hn-l)a"], (24) 



2 

and the forward and backward difference operators are, respectively, 

Af{x) = fix + 1) - fix), V/(x) = fix) - fix - 1). 

It can be proved [8] that the differences of Vnix), namely A"^y„(x) = v^nix) satisfy 
similar equation of hypergeometric type: 

aix) AV V^nix) + Tmix) AVmnix) + l^mnVmnix) = (25) 

with 

Tmi^) = ^i^ + m) + aix + m) — aix), 
l^mn = A„ - = -in - m) It' + ^ ^ — ^cr" ) , m = 0, 1, . . . n - 1. 



We call the polynomials Vmnix) the generalized classical orthogonal polynomials of 
discrete variable, among them we find the Hahn, Chebyshev, Meixner, Kravchuk and 
Charlier polynomials. 

The classical orthogonal polynomials of discrete variable satisfy an orthogonality prop- 
erty with respect to the weight function p(x) 

6-1 

^ yeix)ynix)pix) = 5endl. (26) 

x=a 

Similarly the generalized classical orthogonal polynomials of discrete variable satisfy 
the orthogonality property 

6-1 

^Vmeix)Vmnix)pmix) ^ 5endl,^, (27) 
x=a 

where and are normalization constants. It can be proved [9] 

/n-l \ -1 /n-l > 



\k=m / \fc=0 
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If we eliminate d„„ in the above equations we get 



m— 1 



d'mn = dl n /^fc- (28) 



k=0 



The generahzed clasical orthogonal polynomials of discrete variable can be calculated 
from the weight function Pm{x) with the formula [9]: 



A R 



vU^) = ^!^V"- {pn{x)} , (29) 



. m-1 , , 

n\ -n- / 



where 

(n-m)! \^ ri + Zcy 

(31) 

The leading coefficients of the classical orthogonal polynomial of discrete variable 
yn{x) = OnX" + hnX^~^ + . . ., are given by [10] 



fc=0 



from which it follows AnnBn = Ti!a„. 



We have now all the necessary ingredients to construct the raising and lowering opera- 
tors for the generalized orthogonal polynomials of discrete variable in analogy with those 
of continuous variable. Prom (29) we have 



Wi(^) = ^"'"^;^"^' V"-+MPn+i(a:)} = 
-^'"'"^'^"+-V'^-'"{r„(x)p„(x)} = 



Prom (29) we have 



Am,n+i |^^(^) v"-™p„(a;) + (n - m)r;V"-'"-Vn(a; - 1)} (33) 



V— ' - 1)} = ^W^A-«,„(. - 1) = = V..„(x) 

Substituting the last two expressions in (33) and using (30) we obtain 

Bn+i ( n+1 Xn+rn "-/n /n ^+1 
^^m,n+lO^) = < —. ^rn{x)Vmn{x) — j —T^a{x)VVmn{x) 

Bn yn + 1 —m n + m An n + l — mA„ 



(34) 
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that raises in one step the order of the generahzed polynomials in terms of the polynomials 

Vmn{x) and VVmn{x) 

In order to evaluate the lowering operator we calculate a recurrence relation for the 
generalized orthogonal polynomials of discrete variable. We write 



n+l 



fc=0 
1 

Ckn = -^^Vmk{x) X Vmn{x)pm{x). (35) 

As in the case of the continuous variable Cfc„ = 0, if A; + 1 < n. Hence 

{x) +$nVmn{x) + ) (36) 

where O;^ C^+l,n) /^n (^nni Tn C^— l,n 

Prom (35) it follows that dl^^Ckn = dl^^Cnk 

Since q;„_i = Cn,„_i, % = c„_i,n, if we put A; = n- 1, we get Cn_i,nC?^,„_i = c„,„_id^_„, 
hence 

~ m,n 

""rTijn— 1 

Introducing the expansion yn{x) = a„a;" + bnX"'~^ + ... in (36), comparing the coeffi- 
cients of the highest terms, and using 

777 

A'"^" = n(n - 1) . . . (n - m + l)x"-'" + -n(n - 1) ... (n - m)x'^-'"+^ + . . . 

we obtain 

a„(n — m + 1) = Q;„a„+i(n + 1) 

TTZ in 
ann{n -m)— + bn{n - m) = Q:„6„+in + /SnOn'T' + a;„a„+i(n + l)?^^- 

From these relations and (32) wc obtain 

an n — m + 1 Bn A„ 2ri 2n + 1 n — m + 1 

(X = — 

a„+i n + l Bn+i n \2n+i n+l 

~ hn {n- m) hn+i n-m + l m 

Pn— --; Woj 

ttn n dn+l n+l 1 



(37) 



Hence (34) can be written 



ttn^ Vm,n+l{.x) = -^^^ ^'^^^K mn{x) - a{x)VVmn{x) (40) 



Inserting (36) in (40) we get 



7n ^ 'ym,n-l(2^) = ""^7^ ""T^ 'ymn(2^) + ~ '"rnn(2;) + Cr{x) V'y^n(x) (41) 



The right side of (40) and (41) can be considered the raising and lowering operators 
with respect to the index n for the generahzed orthogonal polynomials of discrete variable 
on homogeneous lattice. 

All the constants an, Pn,Jn, ^n,Tn can be calculated from the tables [11]. Obviously, 
when m — 0, din, Pn,! become, respectively, an, (3n,ln- 

Now we define the orthonormal function of discrete variable 



Using the identity ^P™-^ J = ^ , and the properties of the backwards 

operator we get 



_| Tm-l(3^) 0mn (^) (42) 

^/a(x) + ^/^~Jx)T~a{x) ^/a{x) + Tm-i{x) 

Multiplying both sides of (40) by Pm{x) and inserting the value d'^^yj PnJyx) 

Vvr„n(2;) obtained in (42), we get 

^ ^J'"'^V m,n+l(a;) = L+ (x, n)0^„ (x) = 

, \n+m Tn{x) , ^/a{x) T^-l(x) | ^ 

^ + ^ < Va(x) + V(^(x)+T^_i(x) J 
- •\/(7(x) - (t„_i(x) + (7(x)) V0^„(x) (43) 

Similarly 

In ^ ^J'"~V m,n-l(a^) = L~(x, n)0^„(x) = 

ZTi (J'm,n 

n + m 2n\ J ./a{x) + ^a{x) + Tm-i{x) J 

Vo-(x) {Tm-l{x) + a{x)) V(j)mn{x) (44) 

The expressions (43) and (44) can be considered the raising and lowering operators, 
respectively, for the generalized orthonormal functions on homogeneous lattice. These 
operators are mutually adjoint with respect to the scalar product of unit weight. 

Notice that in (43) and (44) the last term is proportional to V0mn(-s), which in the 
continuous limit becomes the derivative ^l^'mni^)- 

As in [5] we can factorize the raising and lowering operators as follows: 

L~{x, n + l)L~^{x, n) = ix{n) + /x(x + 1, n)H{x, n) 
L'^{x, n)L~{x, n+ 1) — ii{n) + ii{x, n — l)H{x, n + 1) 

where 
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x,n)^ a{x) 

and H{x, n) is the difference operator derived from the left hand side of (25) after substi- 
tuting 4>rnn{x) instead of Vmn{^)- 



4 Raising and lowering operators for classical ortho- 
gonal polynomials of a discrete variable on nonuni- 
form lattice 

Let y{s) a function of discrete variable satisfying the difference equation with respect to 
the lattice function x{s) 

where a{s) = a [a:(s)] , r(s) = r [^(s)] are functions of x{s) of at most of second and 
first degree, respectively. 

It can be proved [12] that the functions t'fe(s) connected with the solutions y{s) by the 
relations 

t^fe(^sO= f ''''11 , v,{s)=y{s) (46) 



Aa;fe_i(s 
k 
2 



Xfe(s) = 1, A; = 0, 1,2... 



satisfy the difference equation 



where 



, , (7(s + A;) - a{s) + r(s + k)^x {s + k - 1/2) 
= A.(s+(k-l)/2) 

''-A+gl^^A + i:.; (49) 

m=0 ' m=0 

provided the lattice functions x{s) have the form 

x{s) = cis^ + C2S + C3 or (50) 
x{s) = ciq' + C2q~' + C3 (51) 

with ci, C2, C3, q, arbitrary constants. 

When //fc = for A; = n in (47) f„ = const. It can be proved that when k < n, Vfe(s) 
is a polynomial in Xk{s) and in particular for A; = 0, Vo{s) — y{s) is a polynomial of degree 
n in x{s) satisfying (45). 
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An explicit expression for A^, when = 0, is given by 

shna; f , , ^, , lsh(n — l)uj ,,1 , ^, 

where a; = |lng, or g = e^'^. For the square lattice (50) lo = 0; and for the q'— lattice (51) 
we have 

shncj g"/^ — ^ 

The polynomials solutions of (45) satisfy the following orthogonality condition with re- 
spect to the weight functions p{s), namely, 

^"^ / 1\ 

ye{s)yn{s)p{s)Ax ( ^ - 2 ) = '^^"^n (53) 

s=a ^ ^ 

Similarly for the differences of the polynomials Vnis), namely, 

Vmn{s)=A^"'^[yn{s)]^Am-lAm-2...Ao[yn{s)], Afe = ^ 



Axk{s) 
it holds 



b-k-l , 

^ Vynl{s)Vran{s)Pra{s)AXyn ( * ~ 3 ) " 



(54) 



m 



where Pm(s) = p(s + m) H cr(s + 

i=i 

It can be proved that the normalization constants satisfy 

(n-l \ ~^ /n-1 

k=m ) \fe=0 

from which can be eliminated: 

m— 1 



^Ixn = n (^^) 
A;=o 

A particular solution of (45) when A = A„ is given by the Rodrigues type formula 

A solution of (47) when /x^ is restricted to A„, namely, Prnn = Pmi^n) = ~ ^m, 
0, 1, ... n — 1, is given by: 

"^"^ ^ P^(^) " ^^"^ ^^ P„.(^) Va;^+i(5) Va;„_i(.)Va;„(.)^^"^'^^ ^''^ 



where 

fc=o 



m—l r Ti m— 1 ■> 

M' TT ^ 

[n — m] ! ^J- [n + k] 
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Formulas (56) and (57) can be written in terms of the mean difference operator [13] 
Sf{s) = / (s + i) - / - i) = A/ (s - i) = V/ (s + i), tliat is to say, 



yn[s) 



Bri 



A R 

Pm{s) 



5x{s) 



(-1) 



6x {s + f ) 



/ n m\ 
(^-2+2) 



(59) 
(60) 



In order to obtain the raising and lowering operators of the classical orthogonal poly- 
nomials on non-homogeneous lattice, we apply the Rodrigues formula (56) 



yn+i{s) 



n+1 



p{s) '^'^^^ ^^"^^^'^^ p{s) Vx,{s) " " " Vx„+i(.) 



{Pn+l(s)} 



Since 



Vpn+l(s) _ Apn+l{s - 1) _ A {a(s)pn(s)} 

Vxn+i{s) Aa;„+i(s-l) Ax„ (s - i) 
using (59) we have 

S 



Tn{s)pnis), 




^^V(") W (s)o (s)\ = 



Sx{s) 



{r„(.-|)p„(.-|)} = 



6 



6x{s) 



Pn [S 



2) 



shnuj , 



Sx {s - i) 



n— 1 



n 1 

Pn\S---- 



(61) 



The last step can be proved by induction for both cases of x{s) on non-homogeneous 
lattice (50) and (51). 

First of all, we transform the properties of these functions [14] given by 



x(s + n) — x(s) 



shnu 
shcj 



Vx ( s 



n + 1 



x{s + n) + x{s) — chncu x{s) + const, 
into the difference relations 



Sx { s — 



shnuj 
shtc; 



6x { s + 



Sx I s 



ox Is + ■ 

- ^Sx (^-5 + ~'~ '^'^ ('^ ~ 2^) } ~ dx{s) 
Suppose it is true that for any two functions of discrete variable it holds 



(62) 
(63) 



6 



6x{s) 



{/(%(«)} = / (5 + ^) ( 



6 



+ 



6x{s) 

shu Sx (s + ^) I 5x (s - I) 



9is) + 

6 



n—l 



9 s 



+ 
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Then using the properties of the mean operator we have: 

Sx{s) [Sx{s-l)) n V 

shnu;Sf{s + ^) 8 ( S V"' / I' 
^ sho; Sx{s + ^)Sx{s) y5x{s-l) J ^ V 2j+--- 

The second and third term on the right hand side can we written: 



Sx(s + ^) [ 5x(s) sho; Sx{s) \ - I) I \ 2 

Using (62) and (63) the expression between curly brackets is equal to sh(n + l)a;/sha;, 
therefore 

Ji+l / I 1 \ / f \ n+l 



5x{s) 



{/(%(.)} = + 9(s)+ 



sh(n+l)^^/(s + f) / 6 \ 

sha; Sx{s + ^i)[6x{s-^J ^\ 

as required. Substituting f{x) = t„ (s — |) and g{s) = p„ (s — |), the terms of lower 
degree become zero, due to the properties of function t„(s). Therefore (61) is proved. 
Using (60) for m = 1 we have 

Vyn{s) Ay^is-l) A^^^Bn ( S V ' f n 1 

- -Vin{s-1) = — 7T^-7 TY Pn{s 



Vx{s) Ax{s-1) ^ \6x{s-l)J ^"V 2 2 

Therefore, (61) can be written 

Alvarez- No darse and Costas-Santos [15] [16] have given the same formula for the lattice 
(51). Here we have proved the similar expression for both cases (50) and (51). 

Prom (64) we can calculate the raising and lowering operators. Instead, we proceed 
to the general case in section 5, and then take the value m — 0. 
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5 Raising and lowering operators for generalized cla- 
ssical orthogonal polynomials of discrete variable 
on non-uniform lattice 

Prom (57) and (60) we obtain 

Vm,n+l{s) = V„+i {Pn+l(s)} = 

Prn V -5 j 

= ^l!!-±l^V(— ) {rn{s)pn{s)} = 

Pm\S ) 

/ \ n—m 

Am,n+lBn+l I ^ \ \ f U - m\ f U - TJl 



Tn \ S ] Pn \ S - 



= p„ (.'-?) IsmJ v-y^-i)"'^'-!)) 

With respect to the new variable s' = s + y, this expression can be easily differentiated 
as in (61) giving 

{/ \ n—m 

5 \ f n-m\ 

(\ n—m— I 
S \ I n m — 1 



From (60) we get 

Vx{s) ~ Ax{s - 1) 



= Vm+l,ni^ - 1) = 



n—m—1 



Am+l.n Bn I 6 \ f U TTl - I 



p^+i(s-l) ^"V^ 2 ' 2 

Using this result and the values for Am,n given in (58) we get the raising operator for 
Vmn{s), namely, 

Vm,n+i{s) - \[n+l-m][n + m]K ^"^"^^^-"^"^^ 

[n + 1] N,;,(,)^I^\ (65) 



[n + l-m\Xn " Vx(s) 

with [n] = corresponding to all values of lattice functions x{s) given in (50) (51). 
In order to construct the lowering operator we use the recurrence relation 

Xm{s)Vmnis) = anV„i,n+lis) + PnVmn{s) + %Vm,n-l{s) (66) 

where Xm{s) = x [s + f) and Vmn{s) = A(™)t/„(s) 

We introduce the expansion yn{s) = anx" {s) + bnx"'~^ {s) + . . . in the recurrence relation 
(66). We have two cases [17] 
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a) Quadratic lattice: x{s) — s{s + 1). 

A("^)x"(s) = n(n - 1) . . . (n - m + + 
+ - 1) . . . (n - m)(2n - 2m + 

which after substitution in the recurrence relation (66) gives 

an n-m + 1 

«n = — (67) 

ttn+l n+l 

p n-m _ bn+i n-m + 1 3_ 

^ an n dn+i n + l 12"* 
_ ctn-i n — m , , 

In — -l2 l^yj 

b) Exponential lattice x{s) = Aq^ + Bq~^ 

A('")x"(s) ^[n][n-l]... [n-m + 1] <""""'(«) + C x'^^'^'^is) + ... 
which after substitution in the recurrence relation (66) gives 

an [n-m + 1] 

«n = -r — -T-^ — ( '0) 

On+i [n + 1\ 

_bn [n - m] _ bn+i [n-m + 1] 

Pn — r -1 Fill V / 

a„ [n] a„+i [n + 1] 

a„_i [n-m] /^..x 



Since a„ = J] (^"[^rrt j ^® obtain 



Bn \n [2n] [2n + l] [n-m + 1] 
Bn+l [n] A2n A2„+l [n + IJ 

which after substituting in (65) gives 

an ..,„H-i(.) = 1^ ^...„(.) - o{s)^^ (74) 
Inserting the recurrence relation (66) in (74) we get 

The last two equation can be considered the raising and lowering operators of gener- 
alized orthogonal polynomials on non-uniform latices for the functions (50) and (51). In 
the first case the parameter [n] should be taken as n. 

In order to complete the picture, we define an orthonormal function 



^mn{s) = dj^s/p^js) Vmn{s) (76) 
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Using the properties of the difference operator and the identity 

Vpm(s) _ Tm-l(s) AXm-1 {s - |) 

Pm{s) a{s) + Tm-l{s) AXm-1 (s - |) 



(77) 



we get 



y cr(s) + r^_i(s) Aa;^_i [s - ^) 
1 

+ , =x 

^/a{sj + yo-(s) + r„_i(s) Aa;^_i (s - I) 



^C7(s) +T^_i(s) Aa;^_i (s - i) 



Multiplying both sides of (74) by <^mJi\/ Pm('S) and subtituting the value d^jj-^ Pm(s) 
^Vmn{s) obtained in (78) we get 



Oin ^ ^[^£!±1 n^,j^i(s) = L+(s,n)Q^n(s) = 



K+m T-„(s) ■v/o-(g) T^-l(g) Va^m (g + |) . , . 



- )j<^isMs) + r_,(.) Ax_, - (79) 
Similarly 

7n ^ "^7"'^ f^m,n-l(s) = -^'"(s, w)f^mn(s) = 

(T(s) Tm-li-s) VXm (s + I) 



+ ^] ^" 2n V ; ^a(.)+w(.) Aa;^(.-i) V^(-) 
nmnis) + y^a(5)a(5) + r^-i(5) Ax^^-i (^s - ^^^^"^^^^^ (80) 

The last two expressions can be considered the raising and lowering operators for the 
generalized orthonormal functions on non- homogeneous lattices of the type (50) and (51). 
It can be proved that these operators are mutually adjoint with respect to the scalar 
product of unit weight. 

As in the previous sections we can factorize the raising and lowering operators as 
follows: 

L~{s, n + l)L"'"(s, n) — ix{n) + u{s + 1, n)H{s, n) 
L'^{s, n)L^{s, n + 1) = /x(n) + u{s, n — l)H{s, n + 1) 
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where 



^2n ^2n+2 



[2n] [2n + 2] 



u{s, n) 




and H{s,n) is the difference operator derived from the left side of (47) after substituting 
^mnis) instead Vmn{s) given in (76). Notice that the expresions for the factorization of 
the raising and lowering operators becomes the same expresions (32) and (33) given in 



We have developped the construction of raising and lowering operators for classical OP 
of discrete variable on non-homogeneous lattice extended also to the generahzed OP on 
homogeneous and non-homogeneous lattice. 

In the last case (generalized OP) the raising and lowering operators can be defined 
with respect to the index n, the order of the OP, or with respect to the index m, the order 
of the difference derivative of the generalized OP, or both. 

In our work we have taken into account only the index n, although we have suggest how 
to complement the calculus with the index m. We have also introduced the orthornomal 
functions of unit weight, more suitable to quantum mechanical applications. 

Our presentation leads to an easier way for the continuous limit (compair with a 
different presentation in [5]). 

We have already worked out some physical application of raising and lowering opera- 
tors on homogeneous and non-homogeneous lattice. For instance, the quantum mechanical 
models for the harmonic oscillator in one dimension (Kravchuk OP), the hydrogen atom 
(generalized Meixner OP) [22], the Heisenberg equation of motion on the lattice (Hahn 
OP) [23] Dirac and Klein-Gordon equation on a homogeneous lattice (discrete exponential 
function). [24] [25] 

Finally the connection between OP on non-homogeneous lattice and the 3nj-Wigner 
coefficients and its application to spin networks models in quantum gravity are now in 
progress. 
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